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We present the results of a neutron scattering study of the high energy phonons in the super- 
conducting graphite intercalation compound CaC6- The study was designed to address hitherto 
unexplored aspects of the lattice dynamics in C&Cq, and in particular any renormalization of the 
out-of-plane and in-plane graphitic phonon modes. We present a detailed comparison between the 
data and the results of density functional theory (DFT). A description is given of the analysis 
methods developed to account for the highly-textured nature of the samples. The DFT calculations 
are shown to provide a good description of the general features of the experimental data. This is 
significant in light of a number of striking disagreements in the literature between other experiments 
and DFT on CaC6- The results presented here demonstrate that the disagreements are not due to 
any large inaccuracies in the calculated phonon frequencies. 

PACS numbers: 71.20.Tx,63.20.kd,74.25.Kc,78.70.Nx 

Keywords: graphite intercalates, superconductivity, electron-phonon coupling, inelastic neutron scattering 



I. INTRODUCTION 



Superconductivity has been known in graphite inter- 
calation compounds (GICs) for over 40 yearsP Until re- 
cently, however, the maximum ambient pressure super- 
conducting transition temperature was limited to about 
2 KP The discovery of superconductivity in CaC6 at the 
enhanced transition temperature T c = 11.5 K changed 
this and reignited the interest in superco ndu ctivity in 
graphite intercalation compounds (GICs)P^ Although 
early reports suggested that the superconducting mech- 
anism involved coupling via acoustic plasmonspthe cur- 
rent consensus is that the superconductivity is close 
to s-wave symmetrjEHSl and phonon mediatedP^ Hav- 
ing said this, some controversy remains as different ex- 
perimental and theoretical techniques implicate differ- 
ent phonons and electrons as being primarily respon- 
sible for the coupling. Experiments have revealed a 
Ca isotope coefficient of 0.5, which is consistent with 
coupling to only the low energy Ca phononsP^l Den- 
sity functional theory (DFT), however, predicts that in- 
plane G&xy vibrations and out-of-plane C 2 phonons con- 
tribute roughly equally, coupling electrons via transi- 
tions involving the interlayer bandP Information about 
electron-phonon coupling can also be inferred from spe- 
cific heat of CaC6 measurements A detailed compari- 
son between the experimental results and calculations^ 
finds that the electron-phonon coupling to the C z modes 



is stronger than DFT predicts. Further adding to the 
inconsistencies, a recent angle-resolved photo-emission 
(ARPES) measurement proposes a entirely different pic- 
ture, as it finds sufficient coupling between the in plane 
C xy modes and the n* electrons to account for T c with- 
out other contributions.^ However, the extraction of this 
coupling is made difficult as the bands are predicted by 
DFT to curve in that region.^ A first step to resolving 
these problems is to measure the phonons to ensure that 
they are accurately predicted by DFT. 

Due to the highly anisotropic structure of CaC6, its 
phonons can be divided into different groups based on 
the predominant atomic motion, which occupy differ- 
ent energy windows. The Ca atoms are much heavier 
and more weakly bonded than the C atoms, so the Ca- 
related phonons occur at the lowest energies typically 
below ~ 40meV. The graphene sheets have strong pla- 
nar bonding, meaning that the out-of-plane C z modes 
fall between 40-80 meV, while the in-plane C xy modes 
are higher energy and dominate the spectrum from 100- 
180 meV. To date the main phonon dispersion measure- 
ments in CaCg were performed using inelastic x-ray and 
neutro n scat tering and concentrate on phonons below ~ 
40meVPEH. Although some phonon modes of the mea- 
sured dispersi on ar e well reproduced by the DFT calcu- 
lations of RefspEH, other modes are substantially under- 
estimated in energy. 7 Furthermore, an anti-crossing of a 
longitudinal acoustic mode with an unknown mode was 
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discovered,^ which was tentatively ascribed to vacancies 
in the Ga.Ce crystal. The higher energy phonons in CaC6 
have only been measured using Raman scatteringPSEU 
Raman scattering normally only measures phonons at 
small momentum transfer near the center of the Bril- 
louin zone. Furthermore, phonons measured by Raman 
scattering are not representative of the phonons across 
the Brillouin zone, as they are subject to non-adiabatic 
effects , whi ch lead to large energy shifts at the zone 
center PH 

For these reasons we have performed a neutron scatter- 
ing measurement of the high energy C xy and C z phonons 
in CaCg. There are two principle aims: to compare ex- 
periment measurement and theoretical predictions for the 
high energy phonons in CaCg and to search for any tem- 
perature dependence of the phonons, which if present, 
might imply strong electron-phonon coupling. We note 
that this is the first neutron scattering measurement of 
the high energy graphitic phonons in any GIC. Section |TT| 
describes the experimental methods. Macroscopic GIC 
samples are typically highly textured, due to the texture 
of the starting graphite, which means that an average of 
the scattered intensity is performed over a large swath 
of reciprocal space when the neutron spectrometer is set 
for nominal energy and wavevector transfers. Section [TTT| 
outlines calculations to account for the sampling of the 
Brillouin zone to allow a detailed comparison between 
DFT and experiment. We discuss the good overall agree- 
ment between experiment and theory and how it could 
be further improved in section [TV} Section [V] outlines the 
implications of our work in light of the current literature 
on CaCg. 



II. EXPERIMENTAL METHODS 

Samples were synthesized using the Li-Ca alloy method 
from ~ 50 strips of highly orientated pyrolytic graphite 
(HOPG)J2^In this method the Li atoms initially diffuse 
into the graphite to form LiCg , opening the graphite lay- 
ers, after which the Ca is able to replace the Li to form 
the final sample. The resulting 700 mg of sample was 
sealed within an aluminium can to protect it from react- 
ing with air, with the strips mounted on an aluminium 
plate with the sample c-axes aligned with one another. 
Thin aluminium foil was used to hold the sample in place, 
to avoid the possibility of glue reacting with the sam- 
ple. The sample was measured using INI at the Institute 
Laue-Langevin (ILL). This instrument is connected to 
the hot source of the ILL reactor, which provides neu- 
trons of sufficiently high energy to access the 200 meV 
energy transfers required in this study. A Cu (2 2 0) 
monochromator was chosen and scans were performed at 
a fixed final energy of 100 meV, and an erbium filter was 
inserted before the double-focusing Cu (2 0) analyzer to 
absorb spurious elastic signalsP^ The sample was loaded 
into a closed-cycle refrigerator to control the sample tem- 
perature. Data were taken at 5 and 20 K and summed, 
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FIG. 1. (Color online) (a) A measurement of magnetization 
M against temperature T displaying the sharp superconduct- 
ing transition in C&Cq. (b) The Brillouin zone of CaC6- (c) 
A (0 L) diffraction scan on the sample measured using INI. 
Tick marks denote the positions of the C&Cq Bragg peaks. 



taking into account changes in the Bose-Einstein factor, 
as these spectra proved to be identical. In this experi- 
ment we use the hexagonal representation of the CaCg 
unit cell with a = b = 4.333 A, c = 13.572 A, a = f3 = 90° 
and 7 = 120°, because this distinguishes the quasi-two 
dimensional nature of the sample more clearly than the 
rhombohedral unit cell. Due to the nature of the starting 
graphite (HOPG), the sample consists of crystallites of 
size ~ 1 /impSl which are aligned in the c-direction and 
powder-like in the ab-plane. Only two components are 
therefore required to describe the scattering vector Q: 
Q r and Q z which are the in-plane and out-of-plane com- 
ponents measured in reciprocal lattice units. 

Figure [ljc) shows a Q z (0 L) elastic diffraction scan 
performed on the sample using INI. All the principal 
peaks are indexed to CaCg denoted by the markers. Only 
very small peaks are seen coming from the aluminium can 
and LiCg and graphite impurities. The graphite and LiCg 
impurities are estimated to each account for ss 5 % of the 
sample by volume, which is consistent with a number of 
published works EEH This is further backed up by Fig. [l] 
(a), which displays a sharp superconducting transition at 
11.5 K obtained from a sample made by the same method. 
A rocking curve performed on the (0 15) Bragg peak 
shows that the c-axis mosaic of the total sample can be 
described by a Gaussian (e~ x / 2(J ) with a = 4.7°. Due 
to the inelastic cross section being inversely proportional 
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to the large energy transfer,^ the neutron flux at the 
detector is relatively small and long counting times (1-2 
days) were required to obtain sufficient statistics for the 
analysis. 

It is instructive to consider the relative neutron scat- 
tering efficiency of the C and Ca atoms within the sam- 
ple, which have scattering cross sections oc — 5.55 and 
oca — 2.78 barns and masses Mc = 12 and Mc a — 40 
atomic mass units. The phonon scattering intensity is 
directly proportional to the cross section and inversely 
proportional to the mass. Consequently when the sam- 
ple stoichiometry is accounted for, the ratio of C to Ca 
phonon scattering is %acM Ca / {^CaMc) ~ 40. Such a 
high ratio implies that the vast majority of the scatter- 
ing will be from the phonons involving predominantly C 
atom motion. 

The selection rule for neutron scattering from phonons 
can be described by a dot product ( Q- e) 2 , where e is the 
polarization of the phonon! 2 -^ Two reciprocal space loca- 
tions were chosen to take data at with high statistics: 
Q r = 6, Q z = was employed to study the C xy phonons 
between 80 and 200 meV and Q r — 2 Q z = 24 was cho- 
sen to study the C z phonons between 20 and 100 meV. A 
preliminary measurement at Q r = Q z — 24 showed rel- 
atively little phonon structure, so this Q was not chosen 
for the final measurement. 



III. THEORETICAL METHODS 

The textured nature of the sample presents a chal- 
lenge for comparison between theory and experiment. 
Detailed calculations were therefore required to facili- 
tate this comparison. The first step is to compute neu- 
tron spectrum for arbitrary Q. Phonon dispersion were 
calculated using DFT in the linear response, as imple- 
mented in the QUANTUM-ESPRESSO cod<P^, us- 
ing the generalized gradient approximation 30 and ultra- 
soft pseudopotentialsP2 The dynamical matrix was cal- 
culated on a 4 x 4 x 4 k-point grid and then Fourier 
interpolated at any point in reciprocal space. Full details 
are given in Refs. H3 and [HO The Fourier interpolated 
dynamical matrix at a given Q vector is then diagonal- 
ized and the neutron structure factor is calculated along 
the lines of Ref. [35] using neutron form factors and a 
Gaussian resolution, compatible with neutron scattering 
experiments. To account for the texture of the sample the 
structure factor is calculated at several points included 
in a toroidal-like volume in reciprocal space, which is 
most straightforwardly described in spherical polar coor- 
dinates (| Q\, 8, (/)). As the sample is powder-like in plane 
this causes all 9 to be sampled equally. The c-axis mosaic, 
as determined by the rocking curve, is the dominant con- 
tribution to the sampling of <f>, which is represented using 
Gaussian sampling with a = 4.7°. The sampling of | Q\ is 
described by another Gaussian with a — 0.1 A" 1 , which 
is consistent with the width of the in-plane diffraction 
peaks. This volume was sampled randomly and about 
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FIG. 2. (Color online) Comparison of the neutron scattering 
measurements and theoretical simulations of the C xy phonons 
in CaC6 measured at Q r = 6, Q z — from 80-200 meV. (a) 
The neutron scattering signal. Black points represent the 
measured spectrum and the solid line is the result of the sim- 
ulation. An exponential background has been subtracted from 
the data, (b) A cut through reciprocal space in the Q z — 
plane, which depicts the regions that contribute most strongly 
to the peak in the phonon intensity simulation in panel (a) 
from 130-150 meV. Black represents the strongest contribu- 
tion and white represents no contribution, (c) A similar cal- 
culation to (b) for the second peak in the data in the energy 
region 150-190 meV. The spectra are the summed result of 
scans performed at 5 and 20 K. 



1000 Q-points were required to achieve good convergence. 
The instrumental resolution function was calculated us- 
ing the RESTRAX software package,^ which performs 
Monte Carlo neutron simulations of the combined effect 
of INI and the sample. This function was convolved with 
the simulated spectra. Each simulation was checked to 
confirm that it was stable to small changes in the mosaic 
parameters. 



IV. RESULTS 

The C xy modes were measured with Q r = 6, Q z — 
and are depicted by the black points in Fig. [2ja), 
which can be compared to the solid line from the the- 
oretical simulation. We first discuss the two main peaks 
which dominate the spectrum at 137 and 172 meV re- 
spectively. These are reasonably well reproduced by our 
first principles calculations. In particular, the peak at 
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FIG. 3. (Color online) (a) The inelastic neu- 
tron scattering signal for the C z modes in CaCg 
measured between 20-100 meV at Q r — 2 Q z = 
24. Black points represent the measured spec- 
trum and the solid line is the result of the simu- 
lation. An exponential background has been re- 
moved from the data to account for the incoher- 
ent scattering from the Li impurities and the alu- 
minium can. (b) A simulated section of reciprocal 
space in the Q z = 24 plane, which depicts the re- 
gions of reciprocal space which contribute most 
strongly to the peak in the phonon intensity sim- 
ulation in panel (a). The spectra are the summed 
result of scans performed at 5 and 20 K. 



137 mcV is reproduced well in the simulation and the 
peak at 172 meV is predicted to be narrower and at 
about 2 meV lower energy. Further calculations were 
performed to analyze the phonons contributing to these 
peaks. Figure [2jb) and[2jc) indicate the relative contri- 
butions to the scattering intensity for the two principle 
phonon peaks in the ranges 130-150 and 150-190 meV 
respectively. The calculations reveal that different re- 
gions of the Brillouin zone provide the majority of the 
phononic signal. The strong variation in scattering ef- 
ficiency from various small areas of the Brillouin zone 
demonstrates why a relatively large number of sampling 
points (~ 1000) were required to adequately represent 
the Brillouin zone. From Fig. ^b) one can see that the 
~ 137 meV peak comes predominantly from Q = (6 0) 
direction and symmetry related points. We note that 
the six-fold symmetry of these phonons arises because 
for the C based modes, the influence of the intercalant 
superlattice is only a weak perturbation of the six-fold 
symmetry of the graphene layers. In this experiment a 
large fraction of the Brillouin zone is sampled in a rather 
uneven way, making the measurement almost, but not 
quite, a 2D density of states measurement. The intensity 
contributing to the peak at 150-190 meV from is shown 
in Fig. [2](c) this intensity comes from a larger region of 
reciprocal space and also displays six-fold symmetry. 

There are also two minor peaks at 120 and 196 meV. 
The peak at 196 meV is we ll abov e the highest calculated 
phonon frequency in CaC6 EElllSI Since our samples have 
roughly 5% of graphite, this peak could be due to the 
in-plane phonon vibrations of pristine graphite sheets, 
which are indeed near this energy.^ The broad peak at 
120 meV is also not completely reproduced in theoretical 
calculations. A peak is present in the theoretical struc- 
ture factor but at substantially lower energy (113 meV) 
and with a lower intensity. It is important to note that 
the lower intensity is due to the highly dispersive na- 
ture of the C xy phonon dispersion in this region. As 
stated before, the in-plane averaged structure factor is 
related to the 2D phonon density of states. In graphite 
the C xy phonons are very flat near the K point at ~ 120 
meV. Consequently this peak could also be due to pris- 
tine graphite. It would be surprising, however, if a 5% 
content of graphite could give a relatively significant sig- 



nal in neutron scattering. 

The C z modes within CaCg were measured at Q r = 
2 and Q z = 24, where the large Q z component acts to 
enhance the scattering power of the z-polarized phonons. 
In this region of the spectrum, the strong incoherent scat- 
tering from the Li impurities and the signal from the alu- 
minium caused a background, which was modeled as two 
exponential functions and subtracted from the data. The 
data has peaks at 54 and 76meV. It can be seen from 
Fig. |3ja) that the simulated spectrum reproduces the 
general shape of the spectrum quite well, however there 
are some differences. Although the position of the peak 
at 76meV is predicted, its relative intensity is underesti- 
mated. Furthermore both the position in energy and the 
relative intensity of the main peak at 54 meV are overes- 
timated. The fact that the 54meV peak is overestimated 
in energy at 56.4 meV is especially interesting as the en- 
ergy of the C z modes are sensitive to the occupation of 
the 7r* band: as the it* band is populated these mode 
softenP* Furthermore, a softening of the Raman active 
C 2 phonon at q ps has been shown, to correlate with 
the increase in T c in BaC 6 , SrC 6 , YbC 6 , to CaC 6 PThis 
mode softens from 60meV in BaCe to 54meV in CaCg- 
The measurement of the energy of this optical phonon in 
this study is thus important and could suggest that the 
amount of charge in the 7r* band is underestimated by 
the DFT model. Further inelastic neutron scattering or 
inelastic x-ray scattering experiments measuring the C z 
phonons could be particularly interesting. Fig. ^b) dis- 
plays the analysis performed for this scan. The majority 
of the neutron scattering intensity from 30-90 meV comes 
from the (0 2 24) direction and the five other symmetry 
related points. 

The temperature dependence of a phonon spectrum 
can often provide information regarding the origin of the 
phonon linewidth. For example, electron-phonon cou- 
pling such as that seen in Ref. \TE\ typically leads to 
temperature-independent linewidths, while anharmonic 
effects usually lead to a reduction in linewidth at low 
temperatures. Studying such a change in linewidth using 
the triple axis mode of INI is rather limited by the flux of 
the instrument. So the scans in Figs. [2] and [3] could only 
be performed at 5 and 20 K, which proved identical. As 
an alternative approach we chose to look for temperature 
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dependent effects in the phonon spectrum by operating 
INI in the Be filter mode. In this mode a block of beryl- 
lium is placed between the sample and the detector, this 
absorbed all neutrons with energy > 5.2 meV. The large 
aperture means that a large proportion of Q is sampled 
and a larger neutron flux is registered at the detector, 
making a complete temperature dependent study practi- 
cal. 

The beryllium filter spectrometer was used to measure 
the C xy phonons at a range of temperatures between 5- 
300 K. Figure |4]ja) displays two measurements at the ex- 
tremes of the temperature range 5 and 300 K, where the 
spectra has been corrected for the small changes due to 
the Bose-Einstein factor in the neutron scattering cross 
section. The error bars shown are representative of the 
error bars on all points. Within the experimental error 
there is no significant change to the spectrum. A much 
higher resolution study using, for example, inelastic x-ray 
scattering might prove a fruitful method to search for any 
more subtle changes in the spectrum. For this measure- 
ment, the sample was oriented such that Q was approx- 
imately parallel to the c-axis, however, due to the large 
aperture significant in-plane Q was also sampled. The 
measured spectrum has strong similarities to the pub- 
lished density of states calculations from Refs. |9] and 
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FIG. 4. (Color online) (a) The beryllium filter measurement 
of the C X y phonons in CaC6 at 5K (•) and 300 K (o), the 
error bars shown are representative of the error bars on all 
points, (b) The theoretical CaC6 phonon density of states 
(DOS) taken from (i) Calandra and MaurP and (ii) Kim ct 
alP^. (c) The theoretical phonon DOS of graphite^ shown 
for comparison with CaC6- 



I3H reproduced in Fig. [4jb) , but the uneven Q sampling 
leads to some phonon modes appearing more strongly 
in the experimental measurement. This means that this 
measurement cannot distinguish between the two CaCg 
phonon calculations in the literature (Refs. [§l and . 
In Fig. jMc) the calculated phonon density of states of 
graphite — is plotted for comparison with CaCg. We see 
that the main peak in graphite is at higher energy than 
the peak in CaC6, because the Ca atoms act to dope the 
graphene sheets in CaCg which increases the C-C bond 
lcngtrPSl and softens the phonons. Other changes to the 
graphene sheets, such as rolling the sheets to make car- 
bon nanotubes, have a much weaker effect on the high en- 
ergy phonons. For example, the phonon density of states 
for single walled carbon nanotubes the 150-200 meV en- 
ergy range is similar to that of graphite.^ 



V. CONCLUSIONS 



We have performed inelastic neutron scattering exper- 
iments to measure the high energy C z and C xy phonons 
in CaCg. Due to the highly textured nature of these 
samples, detailed calculations were required to make the 
comparison between experiment and density functional 
theory (DFT) predictions. Despite the differences dis- 
cussed, overall for both the C xy and the C z phonons we 
find reasonable agreement between experiment and the- 
ory. This shows that DFT predicts the phonon struc- 
ture of C&Cq reasonably well, which is significant in light 
of several conflicting papers^ 12 * 15 !. The measurements 
presented here imply that experimental and theoretical 
studies need to readdress the electronic structure and 
electron-phonon interaction to definitively explain super- 
conductivity in CaCg. 
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